Connections of geometric measure of entanglement of pure symmetric states to 

quantum state estimation 
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We study the geometric measure of entanglement (GM) of pure symmetric states related to rank- 
one positive-operator-valued measures (POVMs) and establish a general connection with quantum 
state estimation theory, especially the maximum likelihood principle. Based on this connection, wc 
provide a method for computing the GM of these states and demonstrate its additivity property 
under certain conditions. In particular, we prove the additivity of the GM of pure symmetric 
multiqubit states whose Majorana points under Majorana representation are distributed within a 
half sphere, including all pure symmetric three-qubit states. We then introduce a family of symmetric 
states that are generated from mutually unbiased bases (MUBs), and derive an analytical formula 
for their GM. These states include Dicke states as special cases, which have already been realized in 
experiments. We also derive the GM of symmetric states generated from symmetric informationally 
complete POVMs (SIC POVMs) and use it to characterize all inequivalent SIC POVMs in three- 
dimensional Hilbert space that are covariant with respect to the Heisenberg-Weyl group. Finally, 
we describe an experimental scheme for creating the symmetric multiqubit states studied in this 
article and a possible scheme for measuring the permanent of the related Gram matrix. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.67.Mn, 03.65. Wj 



I. INTRODUCTION 

Quantum entanglement is a crucial resource for quan- 
tum computation [l|, Q and other information processing 
tasks, such as quantum tclcportation [31 , superdense cod- 
ing [Jj and quantum key distribution [5| . In the past few 
decades, there have been tremendous efforts in under- 
standing various aspects of entanglement in both bipar- 
tite and multipartite settings. One of the central issues 
in entanglement theory is the characterization and quan- 
tification of multipartite entanglement [f| 0] ■ Among the 
many approaches to the investigation of entanglement, 
several geometrically motivated measures have been pro- 
posed and proved to be useful, such as relative entropy 
of entanglement Q , geometric measure of entangle- 
ment (GM) [10l - [l2| and logarithmic global robustness 
[131 , Il4j . Remarkably, these three measures turn out to 
be related fjj-lll. 



Among the three measures, GM seems to be the eas- 
iest to handle and has thus started to attract attention 
in recent years. Moreover, its applications have gone be- 
yond entanglement theory. GM is closely related to the 
construction of optimal entanglement witnesses ll| , ex- 
perimental estimation of entanglement plj , and discrimi- 
nation of quantum states under local operations and clas- 
sical communications (LOCC) [III EH EH • Recently, GM 



has also been utilized to determine the universality of re- 
source states for one-way quantum computation (23l . |24| , 
and to study generic multipartite pure states as compu- 
tational resources [2(], [25| . In the context of condensed 
matter physics, GM has been demonstrated to be useful 
for studying quantum many-body systems, such as char- 
acterizing grou nd state properties and detecting phase 
transitions |26l - (32j . 

Given the above applications of GM, it is thus inter- 
esting to explore its connection with other important re- 
search fields, such as quantum state estimation [33l |3~4|. 
Quantum state estimation is a procedure of inferring the 
state of a quantum system by general measurements — 
positive-operator- valued measures (POVMs). It is a 
central issue of quantum mechanics and a cornerstone 
of various quantum-information processing tasks, such 
as quantum computation, quantum communication and 
quantum cryptography. In this context, many measure- 
ment and reconstruction methods have been proposed 
to estimate the target state. Among the reconstruc- 
tion methods, an efficient one is the maximum likelihood 
method, which has been widely used in experiments [33 - 

Concerning measurement schemes, mutually unbiased 
bases (MUBs) (4l| - |45j and symmetric informationally 
complete (SIC) POVMs pa - (53 |. which stand for very 
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efficient von Neumann measurements and POVMs [54|, 
respectively, are two focuses in the current research com- 
munity. Great efforts have been directed to solving their 
existen ce p roblem and understanding their structure 
0, IH-El since they are closely related to the 

physics in finite-dimensional Hilbert spaces [45|, [49|, [53| . 
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In this article, we aim to establish a general connection 
between GM and quantum state estimation theory, and 
MUBs, SIC POVMs in particular. The cross fertilizing of 
these geometric ideas may bring insights to both research 
fields. We begin by studying the GM of pure symmet- 
ric states related to rank-one POVMs, and establishing 
a general connection with quantum state estimation the- 
ory, especially the maximum likelihood principle [33l - l37j . 
Based on this connection, we provide a method for com- 
puting the GM of these states and demonstrate its ad- 
ditivity property under certain conditions. In particular, 
we prove the additivity of the GM of pure symmetric 
multiqubit states whose Majorana points under Majo- 
rana representation p55l - [59j are distributed within a half 
sphere, including all pure symmetric three-qubit states. 
We then introduce a family of symmetric states that are 
generated from MUBs, and derive an analytical formula 
for their GM. These states reduce to Dicke states in spe- 
cial cases, which are useful for quantum communication 
and have been realized in experiments (i() 7i2l . Next, 
we compute the GM of symmetric states generated from 
SIC POVMs. This result is then used to characterize 
all inequivalent SIC POVMs in three-dimensional Hilbcrt 
space that are covariant with respect to the Hciscnberg- 
Weyl group [4(H48l |52| |. Finally, we propose an exper- 
imental scheme for creating the symmetric multiqubit 
states studied in this article. 

The rest of the article is organized as follows. In 
Sec. lU we establish a general connection between the 
GM of pure symmetric states and the maximum likeli- 
hood principle in quantum state estimation theory. In 
Sec. IIII1 we prove the additivity of the GM of pure sym- 
metric multiqubit states whose Majorana points are dis- 
tributed within a half sphere, including all pure symmet- 
ric three-qubit states. In Sec. IIV1 we introduce a family 
of symmetric states generated from MUBs, and derive 
an analytical formula for their GM. In Sec. [V] we derive 
the GM of symmetric states generated from SIC POVMs 
and use it to characterize inequivalent SIC POVMs in 
three-dimensional Hilbert space. We discuss experimen- 
tal methods for realizing the symmetrized states of this 
article in Sec. I VII We conclude with a summary. 



II. GM OF PURE SYMMETRIC STATES: A 
CONNECTION WITH QUANTUM STATE 
ESTIMATION 

In this section, we study the GM of pure symmetric 
states related to rank-one POVMs, and establish a gen- 
eral connection with quantum state estimation theory, 
especially the maximum likelihood principle. Based on 
this connection, the GM of many pure symmetric states 
can be computed analytically, and its additivity property 
be demonstrated. In the next section, this connection 
will be used to show the additivity of the GM of pure 
symmetric multiqubit states whose Majorana points are 
distributed within a half sphere, including all pure sym- 



metric three-qubit states. Additional examples related to 
MUBs and SIC POVMs will be presented in the following 
sections to illustrate this general idea. 

A. Preliminary 

The geometric measure of entanglement measures the 
maximum overlap between a given state and the set of 
separable states, or equivalently, the set of pure product 
states, and is defined as [Til EH 

A 2 (p) := max tr(per) = max ($|p|$), (1) 
ctgsep v |*)ePR.o 

G{p) := -2 log AOs). (2) 

Here "SEP" denotes the set of separable states, and 
"PRO" the set of pure product states that fully factorize; 
"log" has base 2 throughout this article. Any pure prod- 
uct state maximizing Eq. ([T]) is a closest product state of 
P- 

For symmetric states, the computation of GM can be 
greatly simplified due to a result in Refs. [20l l63l - [65j : 

Proposition 1 The closest product state to any pure or 
mixed symmetric N -partite state p can be chosen to be 
symmetric; it is necessarily symmetric if N > 3: 

A 2 (p) = max f(g)(^|V((g)|^)) 

|¥>i>....,Ivjv) \jZi / ' 

= max^®"^)®". (3) 

m 

Recently, the GM of pure symmetric three-qubit states 
have been derived based on this observation [jjf|, and a 
class of maximally entangled three-qubit states has also 
been obtained |66l. lo7j . The problem is still open for more 
general situations, although progress has been made with 
respect to the connection between the singular values of 
a hypermatrix and the GM (68|. Below, we shall pro- 
vide many examples where analytical solutions can be 
obtained. 

Another key ingredient in our investigation is the max- 
imum likelihood (ML) principle of quantum state estima- 
tion [33U371 ] . Consider state estimation using a rank-one 
POVM composed of M outcomes that are represented by 
subnormalized pure projectors IL, = \ipj)(ipj\ such that 
TLj = I. If we are given N copies of an unknown input 
state and perform N measurements independently, then 
the outcome statistics obey a multinomial distribution. 
Suppose outcome j occurs rij times for j = 1,2 ... ,M 
(Y]j Tij = N); then the frequency of obtaining outcome 
j is fj = rij/N. In the standard state reconstruction, 
the estimator is obtained by solving the following set of 
equations: 

<tl>iWi) = fi Vi- (4) 

However, such a solution does not always exist. The ML 
principle consists in choosing a state pml that maximizes 
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the likelihood functional C(p) as an estimator of the true 
state [H-[33, 

M 

C(p) = l[p?, (5) 

i=i 

where Pj = (ipj\p\ipj) is the probability of obtaining out- 
come j given the input state p. If there exists a state p s 
that satisfies Eq. that is, the probabilities Pj derived 
from this state coincide with the frequencies fj , then p s is 
also the maximum point of the likelihood functional. In 
general, there is an efficient iterative algorithm for finding 
the ML estimator if the POVM is informationally com- 
plete (IC) A POVM is IC if we can reconstruct 
any input state according to the measurement statistics. 
If the POVM is not IC, the maximum of the likelihood 
functional can still be computed efficiently, but the ML 
estimators are generally not unique. In addition, the ML 
principle is also applicable when the (\ipj)(Tpj |)'s form an 
incomplete POVM, that is J2j IV^K^/I = II < J 
(be sure to distinguish "complete" and "informationally 
complete" ) . 

B. Connection 

We are now ready to show the connection between the 
GM of pure symmetric states and quantum state esti- 
mation theory. Following the above notation, assume 
IV'j'KV'jl = n < /, and the largest eigenvalue of II 
is 1. Define ^{{nj})) as the symmetrized state of the 
product state ®f =1 

M 

|*({n J }))=cP sym (g)|V^)^, (6) 
j=i 

where P sym is the projector onto the symmetric subspace, 
and c is a normalization constant, which can be assumed 
to be positive without loss of generality. The effect of 
the projector P sym is determined by its action on pure 
product states. The action of P sy m on the tensor product 
of N single-particle kets \a,k) is given by 

N N 

P sym (g)|a fc ) = — (£)Kw>> (7) 

fc=i ' <?es N fc=i 

where Sn is the symmetry group of N letters. Now sup- 
pose \ak) is the fc-th member of the multiset consisting 
of m copies of n 2 copies of \tp2) and so on. Define 
A as the Gram matrix of the kets |afc), i.e., 

A jk = ( aj \a k ), j,k = 1,2,..., N. (8) 

The dependence of A on \ipj) and rij has been suppressed 
to simplify the notation. The normalization constant c 
is thus given by c = y/ N\ /pcrm(A) (69| , where perm( A) 



denotes the permanent of the matrix A (for other con- 
nections to permanent, see Ref. HH). 

According to Proposition!]] the GM of ^({n-j})) reads 

A 2 (l*({^}))) = maxl^l®^*^-}))! 2 

\<p) 

/VI M 

= p^T,rn |( ^ >r ' 

N\ 

- ^(A) n T c{p) 

< Nl TT -i 

< -— max p. J 

perm(A) o< Pi <i, 52 jPj =i jj^ J 

An M 

perm(A) JUL h ' ^ > 

Note that, apart from a constant factor, the functional 
to be maximized in the second line of the above equa- 
tion is exactly the likelihood functional C(p) associated 
with the POVM {LL,} if we take fj as "frequency" [see 
Eq. J5J], except that p is restricted to pure states here; 
this is the reason for the first inequality in Eq. ©. The 
"frequencies" fj are called compatible if there exists a 
normalized pure state \tp s ) such that | (y? s IV'j) | 2 = fj Yj; 
that is, the fj's can coincide with the probabilities de- 
rived from a pure state. Note that such a pure state 
is necessarily contained in the cigenspacc to the largest 
eigenvalue of II. The maximum in the last line of Eq. © 
can be obtained if and only if the /^'s are compatible. In 
that case, the likelihood functional C(p) is maximized at 
the pure state \<p s ), and \ip s )® N is a closest product state 

to |*(K})>. 

Moreover, if the maximum of C(p) can be obtained 
at a pure state, which is true if the /-,'s are compatible, 
then the GM of ^({jij})) and any pure symmetric TV- 
partite state |$) is additive; that is, G{\^{\nA))®\^)) = 
G(\V({nj}))) + (see Ref. for a detailed dis- 

cussion on the additivity property of GM). We prove 
this statement in Appendix A. In that case, the GM of 
^({rij})) is equal to its asymptotic regularized quantity 
and gives a lower bound for asymptotic relative entropy 
of entan glem ent and asymptotic logarithmic global ro- 
bustness [201 ] . For convenience, we summarize the above 
observations as the following theorem. 

Theorem 2 Suppose there are M kets \ipj) (not neces- 
sarily normalized) that satisfy Sjli IV'j) (V'j I = H an d 
the largest eigenvalue ofH is 1; ni,n2,... ,Um a-re non- 
negative integers with sum N = JV rij and fj = rij/N. 
Define A as the Gram matrix of the multiset consisting 
of ni copies of , ri2 copies of \ip2) an d so on; and 
M{nj})) = N /^Vperm(A)P sym ®f =1 |^)®"' - Then 
the GM of \^>{{nj})) is lower bounded by 

( N\ M \ 
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the 



bound is saturated if and only if the fj 's are compati- 
ble, that is, there exists a normalized pure state \ip s ) such 
that 



(ii) 



Moreover, if the maximum of C(p) [see Eq. f5|) and 
can be obtained at a pure state, which is true if the fj 's 
are compatible, then the GM of \^({nj})) and any pure 
symmetric N -partite state \<5>) is additive; that is, 

® |<f)) = G(\9({ nj }))) + (12) 

No matter whether Eq. (fTTj) can be satisfied or not, 
the maximum of the likelihood functional C(p) in Eq. Q 
can be computed efficiently with an iterative algorithm 
^ g ives an upper bound for A 2 (|*({n,}))) and 
thus a lower bound for G(\^({nj}))). On the other hand, 
when Eq. (fTTj) is not satisfied, it is still possible that the 
maximum of the likelihood functional is obtainable at a 
pure state. Then the iV-fold tensor product of the pure 
state is a closest product state to ^({n,})), and the GM 
of ^({fij})) and any pure symmetric iV-partite state is 
additive (see Sec. IIIII for examples). In other words, the 
compatibility condition on the fj's is sufficient but not 
necessary for Eq. ([12")) to hold. It is really remarkable that 
we can derive the GM of these symmetric states and its 
additivity property from the property of the likelihood 
functional. 

Besides, the condition on the largest eigenvalue of II 
in Theorem [3] is not essential; it is adopted mainly for a 
closer connection with state estimation theory. If instead 
the largest eigenvalue of II is g > 0, then Theorem [2] is 
still applicable as long as fj is replaced by gfj in Eqs. (fTO)) 
and pT|) . One advantage of this alternative convention 
is that \tpj) can now be chosen to be normalized. A sim- 
ple example in the case of qubits is as follows. Suppose 
there are M normalized kets \ipj) whose Bloch vectors 
are distributed on a circular cone around the z axis, such 
that YljLi IVvXV'jl = M(I + ra z )/2 with < r < 1; and 
|*) = ^/M!/perm(A)P sym 0^ 1 |^->. Then has a 
unique closest product state |0)® M and its GM is 



G(|tt» = - log 



Ml /1 + 7'N- 

pcmi(A) \ 2 



(13) 



III. ADDITIVITY OF GM OF PURE 
SYMMETRIC THREE-QUBIT STATES 

Every pure symmetric A^-qubit state \^f) can be writ- 
ten in the form in Eq. ([B]), that is, 



I*) oc R 



M 

sym (g) 

3=1 



0rij 



(14) 



the Majorana representation [55(. Under this representa- 
tion, each pure symmetric A^-qubit state corresponds to 
A^ points on the Bloch sphere. Following Refs. [58l. [59j. 
the \ipj}'s are called Majorana points of |\E'). Recently, 
Majorana representation has found many applications in 
the study of multipartite entanglement, such as classi- 
fication of entanglement of pure symmetric multiqubit 
states under stochastic LOCC (SLOCC) [56], investiga- 
tion of the GM of these states and maximally entangled 
states among them j57l - l59| . 

As an important application of the theory developed 
in the previous section and the Majorana representation, 
in this section, we prove the additivity of the GM of pure 
symmetric multiqubit states whose Majorana points are 
distributed within a half sphere, which is true for all pure 
symmetric three-qubit states. As we shall see shortly, 
the structure of the state space and its boundary plays a 
crucial role in proving this additivity property. 

Theorem 3 Suppose |\&) is a pure symmetric N-qubit 
state whose Majorana points are distributed within a half 
sphere under the Majorana representation; then the GM 
of and any pure symmetric N -partite state is addi- 
tive. 

According to Theorcm[2l to prove this theorem, it suffices 
to show that the maximum of the following functional can 
be obtained at a pure state, 



A I 



c( P ) = i[((^\p\^)r. 



(15) 



Without loss of generality, we can assume that the Majo- 
rana points of |\&) lie within a half sphere with z > 0, and 
the Bloch vector of p is r = (x, y, z) with x 2 + y 2 + z 2 < 1. 
Then it is straightforward to verify that C(p) is nonde- 
creasing with z, and thus its maximum can be obtained 
at the boundary of the Bloch sphere, that is, at a pure 
state. This completes the proof of Theorem [3] 

For any pure symmetric three-qubit state, the three 
Ma jorana points lie within some half sphere; hence, The- 
orem[3]is applicable. The same is true for any pure sym- 
metric multiqubit state that has at most three distinct 
Majorana points. 

Corollary 4 The GM of any pure symmetric N-qubit 
state that has at most three distinct Majorana points and 
any pure symmetric N -partite state is additive. In par- 
ticular, the GM of any pure symmetric three-qubit state 
and any pure symmetric tripartite state is additive . 

Corollary 5 Suppose |\&) is a pure symmetric N -partite 
state which can be written in the form 



I*) oc R 



sym ' 



(16) 



and this representation is unique up to permutations of 
the IV'jTs and some phase factors. This is well-known as 



where n\ + ri2 + n% = N; then the GM of \^>) and any 
pure symmetric N -partite state is additive. 
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Without loss of generality, we can assume the three states 
\tpj) belong to a three-dimensional Hilbert space. As in 
the case of qubit, they can be seen as three extremal 
points of the eight-dimensional state space whose origin 
is the completely mixed state; however, the boundary of 
the state space is no longer a sphere, and the states on 
the boundary are not necessarily pure. In addition, we 
can find a suitable hypcrplanc passing through the origin 
such that the three points are on the same side of the hy- 
perplane (or on the hyperplane). According to a similar 
reasoning that leads to Theorem [3] and Corollary [H the 
maximum of the following functional 



c(p) = Y[((^\ P \^)y 



(17) 



can be obtained at a state pml on the boundary of the 
state space, whose rank is at most two. If pml is pure, 
then we are done. Otherwise, when p is restricted to the 
support of pml, we have 



£(p)=n(^>i^>p. 



(18) 



where \tpj) is the projection of \ipj) onto the support of 
puh- Now applying the same reasoning that leads to 
Corollary 2] shows that the maximum of C(p) can be ob- 
tained at a pure state. Therefore, the Corollary follows 
from Theorem [2] 

Theorem |3] and Corollaries 0] and [5] provide a method 
for computing the asymptotic GM, which in turn pro- 
vides a lower bound for the asymptotic relative entropy 
of entanglement and the asymptotic logarithmic global 
robustness [2(| • They are also useful in the study of mul- 
tipartite pure states as computational resources, since 
GM and its additivity property are closely related to 
whether these states are universal for quantum computa- 
tion [20|, [23 25]. Corollary @] also implies the multiplica- 
tivity of the output purity of the quantum channels asso- 
ciated with pure symmetric three-qubit states according 
to the Werner- Holevo recipe [2(| [7(| ■ 



IV. SYMMETRIC STATES GENERATED 
FROM MUTUALLY UNBIASED BASES 

To illustrate the general idea presented in Sec. UH in 
this section we consider the situation where the POVM 
can be decomposed into a set of von Neumann measure- 
ments, in particular the scenario where the bases of the 
von Neumann measurements are mutually unbiased [IBj . 
We first consider pure symmetric states generated from 
two bases of the qubit Hilbert space, which reduce to 
Dicke states as special cases. We then generalize the 
idea to higher-dimensional Hilbert spaces and point out 
the role played by mutually unbiascdness. 



A. Generalization of Dicke states 

Given r = (0,0,1) and ri = (sin0, 0, cos 0), let |0), 
|1) denote the eigenbasis of Tq ■ cr = a z , and \9 + ) = 
cos(0/2)|O)+sin(0/2)|l), |0_) = sin(0/2)|O) -cos(0/2)|l) 
the eigenbasis of r% ■ cr. Given four nonncgative integers 
{rijk} = {n o, n i;n 10 , n n }, let N = n 0Q + n i + n w + 
nix, Nj = rijo + rijt, f ]k = njk/Nj (assuming Nj ^ 0). 
Define A as the Gram matrix of the multiset consisting of 
n o,n i,ni ,nii copies of |0), |1), |0+), |0_), respectively; 
and 



|*(0,K fc })) : = 



1 N\ ps 
perm(A) sym 



(JO)®" 00 <8 II)®" 01 
(19) 



Here the permanent perm(A) can be computed efficiently, 
see Appendix B. Note that |^(0, {rijk})} can be seen as a 
generalization of the Dicke state; it reduces to the Dicke 
state when = 0,tt or N = or Ni = 0. The GM 
of the Dicke state has been derived in Ref. [ll|, and its 
additivity property has been demonstrated in Ref. poj . 
When No, N\ 7^ 0, we have the following theorem. 

Theorem 6 The GM of |*(0, {n jk })) is lower bounded 
by 



G(|*(0,{ %fc })»>-lo£ 



Nl 



perm(^4) 



f «oo f "oi fnio fin \ . 

Joo hn J 10 -In J- 
(20) 



the bound is saturated if and only if there exists a qubit 
state I if) such that 



|MO)| 2 = / 00 , |(^|l>| 2 = /oi, 
IW + }| 2 = /io, |(^|0-)| 2 -/ii. 



(21) 



When ^ 0,7r, this condition is equivalent to the follow- 
ing one: 



\h s + fusil < 1, 



(22) 



where hj = fjo — fji, and Sq = (—cot 0,0,1), Si = 
(esc 0,0,0) is the dual basis ofro,ri in the x-z plane. 

The GM of | V I'(0, {njk})) and any pure symmetric N- 
partite state is additive, irrespective whether the condi- 
tion Eq. H21\) is satisfied or not. 



We can briefly show Theorem [6] as follows. First 
Eqs. ([2TJ)) and (|2"Tj) can be derived according to a similar 
reasoning that leads to Theorem [2] Equation (|22|) has a 
nice geometric interpretation under the Bloch sphere rep- 
resentation. Suppose the Bloch vector of \<p) is s. The 
first two equations in Eq. ([2~Tj) restrict s to a plane satisfy- 
ing s ■ Tq = hg, which is perpendicular to ro, and the last 
two restrict s to a plane satisfying s ■ r\ = hi, which is 
perpendicular to n. There exists a pure state satisfying 
Eq. (|2"Tj) if and only if the intersection of the two planes 
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passes through the Bloch sphere; that is, the intersection 
of the two planes and the x-z plane lies within the unit 
circle centered at the origin on the x-z plane. This is ex- 
actly what Eq. (|2"2"|) means. In the special case 9 = §, ^r, 
where the two bases |0), |1) and \9±) = |±) are mutually 
unbiased, Eq. (|22|) simplifies to h% + hf < 1. The addi- 
tivity property of the GM of \^(9, {njk})) follows from 
Theorem [3J This completes the proof of Theorem |BJ 

When noi = n m and Tin = n%o, Eq. (|2ip is trivial to 
satisfy; and the iV-fold tensor product of the two eigen- 
states of <7 y , respectively, are exactly two closest product 
states to (6 , {njk})). Besides, Eq. (f2T)]) reduces to 



G(|tt(0,{n ifc })» 



log 



2 iV pcrm(A) 



(23) 



Remarkably, the GM is completely determined by the 
number of qubits N and perm(A). An interesting exam- 
ple is the state {1, 1; 1, 1})) obtained when noo = 
tiqi = 7)io = = 1- It is a balanced four-qubit Dicke 
state when 9 = 0, 7r, and is equivalent to four-qubit GHZ 
state under a suitable local unitary transformation when 
9 = f,^, according to Refs. HflM El- Calculation 
shows that perm(^4) = [7 + cos(20)]/2, hence, 



In general, it is not easy to determine whether such a 
state exists or not; the structure of the state space plays 
a crucial role here. Here we are content to point out 
a connection with MUBs. Suppose the b bases \e J k ) are 
mutually unbiased, and njk = Tij,fjk = l/dVj,k. Then 
satisfying the set of constraints in Eq. (|25|) amounts to 
the existence of a pure state that is mutually unbiased 
with all states \e k ). 

Theorem 7 Suppose \ei) for j = 0, 1, . . . , 6—1(6 < d+1) 
and k = 0, 1, . . .,d - 1 are b MUBs, i.e., \{e? k \e l m )\ 2 = 
4(1 — Sjj) + 5jj8k t m- Then the N -partite (N = 
dY^j=i n j> n j > 1) symmetric state \*S?(d,{nj})) := 

y/N\/perm(A)P Bym (<g) b j=1 ®^ =1 \4) 9nj ) has GM lower 
bounded by 

G(Md,{n jm > -lo g ^^ y (26) 

the bound is saturated if and only if there exists a pure 
state that is mutually unbiased to all states \e 3 k ). If such 
a state exists, then the GM of \^(d, {nj})} and any pure 
symmetric N -partite state |$) is additive; that is, 



G(|*(0,{1,1;1,1})» = log 



7 + cos(26») 



G(|*(d,K-}))®|$» = GQ*(d,{ nj }))) + GQ*))- 



(24) 



(27) 



The maximum is obtained at 6 = 0,tt, where the two For sample, the inequality in Eqs. ® is saturated 



bases coincide, and the minimum at 9 = f^, where 
the two bases are mutually unbiased. 



B. A connection with MUBs 

To generalize the above idea to higher dimension, let 
|4) for j = 0, 1, . . .,6- 1, k = 0, 1, . . .,d- 1 be bd single- 
particle states such that the d states for given j form 
an orthonormal basis [45[. Let njk be a b x d matrix 
composed of nonnegative integers, N = J2j k n jk> Nj — 
^2 k njk, and fjk = njk/Nj (assuming Nj ^= 0). Define 



2 and the le^Vs 



and (|2T|) is applicable when b 
the eigenbases of Z and X, respectively, where Z and 
X are phase operator and cyclic shift operator, respec- 
tively. They are defined according to their action on the 
computational basis lefe) 



-kh 



Z\e k ) = oo k \e k ), Lo = e 2 ^' d 

X\ek) = |e(fe+l) modd), 



(28) 



Minjk})) : 



1 M p 
perm(A) sym 



i|e£)®"ifc. 



This state reduces to a generalized Dicke state when 
there is only one basis, that is b — 1. According to the 
same reasoning as before, the GM of ^({Hjfc})) is lower 
bounded by 



where we use "mod" to denote the modulo operation. 
Since the eigenbases of X, Z and XZ are mutually unbi- 
ased [HI , the iV-fold tensor product of any eigenstate of 
XZ is a closest product state to \^(d,{nj})). 

If d is a prime power, there exists a complete set of d+1 
MUBs HlHi]. The inequality in Eqs. (gSJ) is saturated 
and (|27p is applicable if |e^.} 's are chosen from 6 bases 
with 1 < b < d from the complete set. However, this 
is not the case if 6 = d + 1, since there is no pure state 
that is mutually unbiased to all states in a complete set 
of MUBs IH. 



G(M{njk}))) > -^g 



( Nl nr 

Vperm(A) 11 J ^ 



SYMMETRIC STATES GENERATED FROM 
SIC POVMS 



the bound is saturated if and only if there exists a pure 
state \ip) such that 



V. 



In a d-dimcnsional Hilbert space, a SIC POVM [46 
|48|. l5ll consists of d 2 outcomes that are subnormal- 



(25) 



ized projectors onto pure states n^- 
j = 1 , . . . , d 2 , such that 



In that case, the GM of [^({n.^})) and any pure sym- 
metric iV-partite state is additive . 



\(^m\ 2 = 



1 + dSjk 



al^-)(^l for 



(29) 
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The condition Y^j Hj = / is already implied by the above 
equation and need not be imposed separately. Most 
known SIC POVMs are generated from a fiducial state 
under the action of the Heisenberg-Weyl (HW) group, 
which is generated by the two operators X, Z defined in 
Eq. ([28]) . A fiducial state of the HW group obeys the 
following equations, 



|<V|x fci z fc2 |v>)| = 



i 



Vd+l 



(30) 



for all 7^ (0,0) mod d. If is a fiducial state, 

then the d 2 states X kl Z k2 \tp) for ki, k 2 = 0, 1, . . . , d - 1 
form a SIC POVM that is covariant with respect to the 
HW group. The Clifford group is the normalizer of the 
HW group that consists of unitary operators. Likewise, 
the extended Clifford group is the larger group that con- 
tains also anti-unitary operators. For any operator U 
in the extended Clifford group, U\ip) is a fiducial state 
whenever is. Fiducial states and SIC POVMs form 
disjoint orbits under the action of the extended Clifford 
group. SIC POVMs on the same orbit of the extended 
Clifford group are equivalent in the sense that they can be 
transformed into each other with unitary or antiunitary 
operations {48l l52j. 

In this section we study the GM of symmetric states 
generated from SIC POVMs, 



l^ IC > 



d 2 \ 



perm(A) 



isym(l^l) ® 



(31) 



Since the completely mixed state is the only state that 
satisfies (ipj\p\ipj) = 1/d for j = l,2,...,d 2 , the GM 
of l^ 10 ) cannot be computed according to Theorem [5] 
In Appendix C, we derive the GM of |4'^ IC ) by virtue 
of the special properties of SIC POVMs; the result is 
summarized in the following theorem. 

Theorem 8 Suppose d 2 normalized states \ipj) for j = 
1,2, ...,d in a d- dimensional Hilbert space satisfy 

d 2 

£.7=1 IV'jXV'jl = dl> an d A is the Gram matrix of the 
\ipj)'s; define |*) = ^(P\/peTm(A)P sym (\ipi) ® • • • ® 
\tpfp.)). Then the following four statements are equiva- 
lent. 

1. The d 2 states \ipj) form a SIC POVM, that is, they 
satisfy Eq. 



2. Up to global phases, the d 2 states \ijjj)® d for j = 
1, ■ • • , d 2 are the only closest product states to |^); 
and 



A 2 m) 



d 2 \ 



(32) 



3. The d 2 states \ip 3 ) m for j = 1, • • • ,d 2 satisfy 

\<M**m\ a = 



(d+ l) d2 - 1 pcrm(A) : 



(33) 



There exists a SIC POVM consisting of d 2 states 
\(fj) for j = 1, 2, . . . , d 2 such that |</?j)® d satisfies 

d 2 l 



(d+ l) d2 - 1 perm( J 4)' 



(34) 



In dimension two, there is only one orbit of fiducial 
states, one of them is given by |46l-|48l| 



14) = v ( 3 + ^)/ 6 i e ") + e ™N ( 3 - ^)/ 6 l ei )'( 35 ) 

Moreover, all SIC POVMs are equivalent to the one thus 
generated. The four states of each SIC POVM form a reg- 
ular tetrahedron when represented on the Bloch sphere. 
Calculation shows that perm(^4) = |, which, together 
with Theorem [51 implies that 

A 2 (|*| IC » = |, G(|* 2 SIC »=log3. (36) 

This result coincides with that obtained by Martin et 
al. [13] and Aulbach et al. [5^]. They also showed by 
numerics that the state generated from the SIC POVM is 
the maximally entangled pure symmetric four-qubit state 
with respect to GM. 

On the other hand, there is a one-parameter family of 
fiducial states in dimension three I46H48I. [521 , 



1 



V2 



(M-^lea)); 



(37) 



and there is a one-to-one correspondence between orbits 
of the extended Clifford group and the parameter t for 
t E [0| §]• There are three kinds of orbits, two excep- 
tional orbits corresponding to the endpoints t = and 
t = ^, respectively, and infinitely many generic orbits 
corresponding to < t < ^ . In addition, there is a one- 
to-one correspondence between inequivalent SIC POVMs 
and the parameter t for < t < f ; the SIC POVMs on 
the three orbits t, =§■ — t, + t, respectively, are equiv- 
alent under unitary transformations. All these inequiva- 
lent SIC POVMs can be classified in terms of geometric 
phases associated with fiducial states (52|. Here we pro- 
vide an alternative characterization in terms of GM. 

Let |^ / f IC (^)) denote the symmetric state related to 
the SIC POVM generated from the fiducial state \tp f 3 (t)) 
under the HW group, A{t) the corresponding Gram ma- 
trix. Surprisingly, the permanent of Ait) can be given 
by a simple formula: perm(A(i)) = |Z [61 — cos(9t)] . Ac- 
cording to Theorem [5] and Eq. @, 



G'l 



\*f c (t))) = log 



16[61 - cos(9t)] 



105 



(38) 



Figure [U shows G(t) := G(\^f G (t))) for < t < f . 
Up to the equivalence class, each SIC POVM is uniquely 
specified by G(|*f IC (f))). There is only a small differ- 
ence among the GM of symmetric states generated from 
different SIC POVMs. The minimum and the maximum 
of the GM are obtained at the exceptional orbits t = 
and t = ?, respectively. Note that the SIC POVMs on 



the orbits t = 
on the orbits t 



3f and t 



are equivalent to the ones 



and t = respectively 
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0.2 0.4 0.6 0.8 1 

t 



FIG. 1: The GM of symmetric states j*f c (t)> for < 
t < j generated from HW covariant SIC POVMs in three- 
dimensional Hilbert space. Each SIC POVM is uniquely spec- 
ified by G(t) up to the equivalence class; recall that there is a 
one-to-one correspondence between inequivalent SIC POVMs 
and the parameter t for < t < ^ 



VI. CREATING SYMMETRIC MULTI-QUBIT 
STATES 

In this section, we describe how to create symmetric 
multi-qubit states mentioned in this article, following the 
approach by Bastin et al. [7ll |. For non-qubit systems, 
a generalization could be made, but we do not know 
whether suitable physical systems exist. 

Imagine we have an atom or trapped ion initially in 
the excited state |e), and it has two stable ground states, 
labeled by |0) and |1), which are connected to the excited 
state via emitting right-circularly (R) and lcft-circularly 
(L) polarized light, respectively. If the detection of the 
emitted photon is R (L) then we know the atom is now 
in the state |0) (|1)). The photon polarization carries 
the which-way information. If we place a polarizer in 
(\R) + \L})/^/2 to erase the which-way information and if 
the detector behind the polarizer clicks, the atom is then 
projected to an equal superposition of |0) and |1). In 
general, if the polarizer is in a\R) + f3\L), after detecting 
the photon, the atom is in the state a|0) + 

Now suppose we have iV excited atoms and N photon 
detectors. A simultaneous detection of N photons will 
project the TV-atom state, depending on the polarizers' 
settings. Suppose the fc-th polarizer is in a.k\R) + Pk\L). 
If the fc-th polarizer and detector sit close to the fc-th 
atom and detect the emitted photon by the nearby atom, 
and if each detector-atom system is far away from each 
other, then after detection of a photon at every detec- 
tor, the iV-atom state is ®k{ak\fyk + /3fc|l)fc). Bastin et 
al. [7l[ proposed the setup such that the N detectors are 
placed in the far-field, roughly at an equal distance to 
all atoms. Because of multipath quantum interference 
(which-path information is erased), after the detection, 
the atoms should be in a state that is invariant under 



permutations, namely, 

|*)ocP sym (g)(a fe |0} fc +/3 fe |l) fe ). (39) 

k 

Thus, using the scheme by Bastin et al. [7l|, the sym- 
metric multi-qubit states discussed in this article can be 
realized in principle, not just mathematical objects. 

For non-qubit systems, one needs to first identify a 
physical system that contains an excited state coupled to 
d longlived sublevels and that the decay into each level is 
associated with a distinguishable which-way information. 
Deleting the which-way information from the decay and 
the which-atom information will enable the creation of 
symmetric multiqudit states considered in this article. 
But we have not yet identified such a system. 

Next, we discuss a possible measurement of the per- 
manent of Aij = (ipi\ipj). Suppose one can prepare the 
following symmetric TV-particle state 

N 

|f)=cP sym (g)|^>, (40) 

3=1 

with the normalization c = \J A r !/pcrm( J 4). If one 
chooses a von Neumann measurement at each particle in 
the basis consisting of \ip) and the remaining orthonor- 
mal states, the simultaneous detection of outcome \tp) at 
all sites has the probability 

VI N 

X ^ 2 = tttIII^)! 2 - ( 41 ) 

perm i4 - LJ - 

3=1 

Since all K^IV^I's are known and \{f) 2 is obtained from 
the measurement statistics, perm(A) can be inferred. 
Ideally, one prefers \(<p) 2 to be as large as possible, as it 
represents the probability of obtaining the desired out- 
come. One can also vary \ip); the probability \(ip) 2 is 
maximized when \<p)® N is the closest product state that 
results in the geometric measure. According to a recent 
work by Martin et al. [57|, the maximal overlap A 2 for 
highly entangled symmetric states decays only inversely 
proportional to the number of qubits, so the proposed 
scheme seems to be feasible even for medium size of sys- 
tems that have been achieved in various physical imple- 
mentations. 



VII. SUMMARY 

We have studied the GM of pure symmetric states re- 
lated to rank-one POVMs and established its connection 
with the maximum likelihood principle in quantum state 
estimation theory. Based on this connection, we pro- 
vided a method for computing the GM of these states 
and demonstrated its additivity property under certain 
conditions. In particular, we proved the additivity of 
the GM of pure symmetric multiqubit states whose Ma- 
jorana points are distributed within a half sphere, in- 
cluding all pure symmetric threc-qubit states. We then 



9 



introduced a family of symmetric states that are gener- 
ated from MUBs and derived an analytical formula for 
their GM. We also derived the GM of symmetric states 
generated from SIC POVMs and used it to character- 
ize all inequivalent HW covariant SIC POVMs in three- 
dimensional Hilbert space. A scheme for creating the 
symmetric multiqubit states studied in this article was 
also proposed. Our studies promise a broad perspective 
of integrating two important research areas in quantum 
information science, namely, entanglement characteriza- 
tion and quantum state estimation. 
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and so on. Then according to Proposition [T] 



A 2 (|*(K}))®|$)) 

= max\{ ( p\® N \(\*({n j }))®\$))\ 

m 



perm(^) \ v ) 
Nl 



N 



)(f\A k ,B k \a k ) Ak )\$)b u ...,b 



perm(^) \ v ) 



k=l 
N 



A k ,B, 



\ a k)A k \ 



fe=l 



\(<P\A k ,B k \a k ), 



\*)B lt ...,B t 



< A 2 (|<&» 



N 



perm(.A) p 



maxJJ ((a k \ Ak PA k \a k )A k ) 



k=l 



A 2 (|1>»- 



max£(/j) 



pcrm(j4) p 
A 2 (|*({n,})))A 2 (|$)). 



(43) 



Here the last equality is due to our assumption that the 
maximum of the likelihood functional C(p) can be ob- 
tained at a pure state. 



Appendix B: perm(A) in Eq. <|19|) 



Note that the entries of the Gram matrix A in Eq. ([19]) 
only take on five different values 0, 1, sin(0/2), ± cos(#/2). 
In contrast with the computation of the permanent of a 
generic matrix, the permanent of A can be computed 
efficiently as follows: 



Appendix A: Additivity of the GM of |*({nj})) 



In this appendix, we prove Eq. (|12[) when the maxi- 
mum of the likelihood functional C{p) can be obtained 
at a pure state. Equation (fT2j) is equivalent to 



A 2 (|*(K})> ® |$)) = A 2 (|v&({n J })))A 2 (|$)). (42) 



perm(A) = \{ n jh \\ £ (-l) 9 " a |c 

\j,k=0 J {a,b,c,f,g} 
Q,2a+2b+2c-2g / n 

x | sin - 

2/ \a, o, noo — a — o 



2> 



raoi 



c, /, noi — c — fj \g,a + c - g, mo — a — c 



a + b — g, f + g — a,n n —b—fj' 



(44) 



According to the definition, the l.h.s is never smaller than 
the r.h.s; so it suffices to show that A 2 (|*({n J }))(K)|$)) < 
A 2 (|*({n i })))A 2 (|$)). Suppose \^({rij})) is shared over 
the parties A\, A2, . . . , Ajy, and |$) over the parties 
B\, £?2, • ■ • , -B/v; suppose \a k ) is the fc-th member of the 
multiset consisting of n\ copies of |^i), n-2 copies of ^2) 



where the summation is restricted to the set of nonneg- 
ative integers {a, b, c, /, g} satisfying the following con- 
straints, 

a + b<n 00 , c + f<n 01 , a + c<n w , b + f<n,n, 
g < a + b, g < a + c, f + g > a. (45) 
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Appendix C: Proof of Theorem [8] 



and [4] can be shown as follows, 



We begin by proving the implication [T}=>[2] According 
to Proposition [T] and the definition of a SIC POVM, 



A 2 (|*» = A 2 (|*i IC ))=max|<^ 

\<p) 



d 2 \ 

7~T\ max 

perm(A) \ v ) 



(46) 



Here the maximization is subjected to the completeness 



d 2 \ 



perm(^4) 



< 



fe=i 

d 2 \ /e£ Wj -k 



2\ d 2 -l 



perm(^4) \ d 2 — 1 

(d + lY^permiA) ' 



condition V - \(M\ 2 = d and the following condition Th ? Equality is saturated if and only if |(^# 



(49) 



k = 



£i<v#i>i 4 



2d 
d+l' 



(47) 



l/(d + 1), Vj, k and j ^ fc; that is, the |^)'s form a SIC 
POVM, which implies [H and |U 

It remains to show the implication [5J=H!I1 



since a SIC POVM is a 2-design [47|: iV'jXV'il ® 
I'^'jXV'jl = [2d/(d+ l)]II sym , where II sym is the projector 
onto the bipartite symmetric subspace. Hence, 



A-(|* 



D) < 



d?\ - 

—-- max pj 

perm(^4) o< Pj <i, £ jPj =d. n 



rf 2 ! 



(d+ l) d -- 1 pcrm(A) 



(48) 



The maximum in the above equation is obtained if all 
Pj's are equal to l/(d + 1) except one of them, which is 
equal to 1. A state can satisfy these conditions if and 
only it belongs to the SIC POVM, hence, the implication 
H=>- M follows. 

The implication[2]=>[3]is obvious. The implication[3]=^T] 



d 2 \ 



perm(A) 
d 2 l 



ni(^i^>i 2 
fc=i 

flfflK^I^ 



perm (^) i = i yii 



d 2 \ 

< j-tt- max 



< 



perm(^) w ^J. 
d 2 \ 

(d+ l) d2 - 1 perm(A)' 



(50) 



Here the second inequality follows from Eq. (|48)). recall 
that the \tpj)'s form a SIC POVM. The inequalities are 
saturated if and only if the \ipjYs are the same as the 
|</?j)'s up to some permutation and phase factors; which 
implies 1. 
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